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Abstract 

This paper deals with numerical solutions to an impulse control problem arising from 
optimal portfolio liquidation with bid-ask spread and market price impact penalizing 
speedy execution trades. The corresponding dynamic programming (DP) equation is a 
quasi- variational inequality (QVI) with solvency constraint satisfied by the value func- 
tion in the sense of constrained viscosity solutions. By taking advantage of the lag 
variable tracking the time interval between trades, we can provide an explicit back- 
ward numerical scheme for the time discretization of the DPQVI. The convergence of 
this discrete-time scheme is shown by viscosity solutions arguments. An optimal quan- 
tization method is used for computing the (conditional) expectations arising in this 
scheme. Numerical results are presented by examining the behaviour of optimal liqui- 
dation strategies, and comparative performance analysis with respect to some bench- 
mark execution strategies. We also illustrate our optimal liquidation algorithm on real 
data, and observe various interesting patterns of order execution strategies. Finally, 
we provide some numerical tests of sensitivity with respect to the bid/ask spread and 
market impact parameters. 
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1 Introduction 



Portfolios managers define "implementation shortfall" as the difference in performance bet- 
ween a theoretical trading strategy and the implemented portfolio. In a theoretical strategy, 
the investor observes price displayed by the market and assumes that trades will actually be 
executed at this price. Implementation shortfall measures the distance between the realized 
transaction price and the pre-trade decision price. Indeed, the investor has to face several 
adverse effects when executing a trading strategy, usually referred to as trading costs. Let 
us describe the three main components of these illiquidity effects: the bid/ask spread, the 
broker's fees and the market impact. The best bid (resp. best ask) price is the best offer 
to buy (resp. to sell) the asset, and the bid/ask spread is the difference (always positive in 
the continuous trading session) between the best ask price and best bid price. The broker's 
fees are the amount paid to the broker for executing the order. The market impact refers 
to the following phenomenon: any buy or sell market order passed by an investor induces 
an adverse market reaction that will penalize quoted price from the investor point of view. 

Market impact is a key factor when executing large orders. A famous worst case example 
is Jerome Kerviel's liquidation portfolio, operated by Societe Generale in 2008. According to 
the report of Commission Bancaire, the liquidative value of Kerviel's portfolio was -2,7G€! 
when the Societe Generale decided to unwind it on January 20, 2008. The liquidation 
was operated during 3 days and led to a supplementary loss of 3,6G€!. Even in regular 
operations, price impact may noticeably affect a trading strategy. On April 29, 2010, 
Reuters agency reports that Citadel Investment Group sold 170M shares of the E*Trade 
stock, and raised about 301M$: this operation led to a price fall of 7,1%. These examples 
explain why measurement and efficient management of market impact is a key issue for 
financial institutions, and the research of low-touch trading strategies has found a great 
interest among academics. 

Most of market places and brokers offer several common tools to reduce market impact. 
We can cite as an example the simple time slicing (we will refer to this example later as the 
uniform strategy): a large order is split up in multiple children orders of the same size, and 
these children orders are sent to the market at regular time intervals. Brokers also propose 
more sophisticated tools as smart order routing (SOR) or volume weighted average price 
(VWAP) based algorithmic strategies. Indeed, one basic observation is that market impact 
can be reduced by splitting up a large order into several children orders. Then the investor 
has to face the following trade-off: if he chooses to trade immediately, he will penalize his 
performance due to market impact; if he trades gradually, he is exposed to price variation 
on the period of the operation. Our goal in this article is to provide a numerical method 
to find optimal schedule and associated quantities for the children orders. 

Recently, there has been considerable interest for this problem in the academic lite- 
rature. The seminal papers [5j and [2] first provided a framework for managing market 
impact in a discrete- time model. The optimality is determined according to a mean- variance 
criterion, and this leads to a static strategy, in the sense that it is independent of the stock 
price. Models of market impact based on stylized order book dynamics were proposed in 
[15], [19] and [9]. There also has been several optimal control approaches to the order 
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execution problem, using a penalizing function to model price impact: the papers [18] 
and [8] assume continuous-time trading, and use an Hamilton- Jacobi-Bellman approach for 
the mean- variance criterion, while [TO], [13], and [IT] consider real trading taking place 
in discrete-time by using an impulse control approach. This last approach combines the 
advantages of realistic modelling of portfolio liquidation and the tractability of continuous- 
time stochastic calculus. In these papers, the optimal liquidation strategies are price- 
dependent in contrast with static strategies. 

In this article, we adopt the model investigated in [TT] , Let us describe the main features 
of this model. The stock price process is assumed to follow a geometrical Brownian motion. 
The price impact is modelled via a nonlinear transaction costs function, that depends 
both on the quantity traded, and on a lag variable 9 tracking the time spent since the 
investor's last trade. This lag variable will penalize rapid execution trades, and ensures 
in particular that trading times are strictly increasing, which is consistent with market 
practice in limit order books. In this context, we consider the problem of an investor 
seeking to unwind an initial position in stock shares over a finite horizon. Risk aversion of 
the investor is modelled through a utility function, and we use an impulse control approach 
for the optimal order execution problem, which consists in maximizing the expected utility 
from terminal liquidation wealth, under a natural economic solvency constraint involving 
the liquidation value of portfolio. The theoretical part of this impulse control problem is 
studied in [11] , and the solution is characterized through dynamic programming by means 
of a quasi-variational inequality (QVI) satisfied by the value function in the (constrained) 
viscosity sense. The aim of this paper is to solve numerically this optimal order execution 
problem. There are actually few papers dealing with a complete numerical treatment of 
impulse control problems, see [6], [H], or [7]. In these papers, the domain has a simple 
shape, typically rectangular, and a finite-difference method is used. In contrast, our domain 
is rather complex due to the solvency constraint naturally imposed by the liquidation value 
under market impact, and we propose a suitable probabilistic numerical method for solving 
the associated impulse control problem. Our main contributions are the following: 

• We provide a numerical scheme for the QVI associated to the impulse control problem 
and prove that this method is monotone, consistent and stable, hence converges to 
the viscosity solution of the QVI. For this purpose, we adapt a proof from [I]. 

• We take advantage of the lag variable 6 to provide an explicit backward scheme 
and then simplify the computation of the solution. This contrasts with the classical 
approach by iterative sequence of optimal stopping problems, see e.g. [6]. 

• We provide the detailed computational probabilistic algorithm with an optimal quan- 
tization method for the approximation of conditional expectations arising in the back- 
ward scheme. 

• We provide several numerical tests and statistics, both on simulated and real data, 
and compare the optimal strategy to a benchmark of two other strategies: the uniform 
strategy and the naive one consisting in the liquidation of all shares in one block at 
the terminal date. We also provide some sensitivity numerical analysis with respect 
to the bid/ask spread and market impact parameters. 
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This paper is organized as follows: Section 2 recalls the problem formulation and main 
properties of the model, in particular the PDE characterization of the impulse control 
problem by means of constrained viscosity solutions to the QVI, as stated in [TT] . Section 3 
is devoted to the time discretization and the proof of convergence of the numerical scheme. 
Section 4 provides the numerical algorithm and numerical methods to solve the DPQVI. 
Section 5 presents the results obtained with our implementation, both on simulated and 
historical data. 

2 Problem formulation 

2.1 The model of portfolio liquidation 

We consider a financial market where an investor has to liquidate an initial position of y 
> shares of risky asset by time T. He faces the following risk/cost tradeoff: if he trades 
rapidly, this results in higher costs due to market impact; if he splits the order into several 
smaller blocks, he is exposed to the risk of price depreciation during the trading horizon. 
We adopt the recent continuous-time framework of [TT], who proposed a modeling where 
trading takes place at discrete random times through an impulse control formulation, and 
with a temporary price impact depending on the time interval between trades, and including 
a bid-ask spread. 

Let us recall the details of the model. We set a probability space (fi, J 7 , P) equipped 
with a filtration F = (J-t)$<t<T supporting a one-dimensional Brownian motion W on a 
finite horizon [0,T], T < oo. We denote by P t the market price of the risky asset, by X t 
the cash holdings, by Yf the number of stock shares held by the investor at time t and by 
@t the time interval between t and the last trade before t. 

Trading strategies. We assume that the investor can only trade at discrete time on 
[0,T]. This is modelled through an impulse control strategy a — (r n Xn)n>i where r\ < 
• • • T n < . . . < T are stopping times representing the trading times and n > 1, are 
T Tn -measurable random variables valued in R and giving the quantity of stocks purchased 
if Cn > or selled if £ n < at these times. A priori, the sequence (r n Xn) m &y be finite 
or infinite. We introduce the lag variable tracking the time interval between trades, which 
evolves according to 

6t = t - T n , T n < t < T n + U 9 Tn+1 =0, U > 0. (2.1) 

The dynamics of the number of stock shares Y is then given by : 

Y s = Y Tn , r n <s< r n+ i, Y Tn+1 = Y Tn + Cn+i, n > 0. (2.2) 

Cost of illiquidity. The market price of the risky asset process follows a geometric 
Brownian motion: 

dP t = P t (bdt + adW t ), (2.3) 

with constant b and a > 0. We do not consider a permanent price impact, i.e. the lasting 
effect of large trade, but focus here on the temporary price impact that penalize the price 
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at which an investor will trade the asset. Suppose now that the investor decides at time 
t to trade the quantity e. If the current market price is and the time lag from the last 
order is 0, then the price he actually get for the order e is: 

Q(e,p,6) = pf(e,0), (2.4) 

where / is a temporary price impact function from R x [0, T] into U {oc}. Actually, in 
the rest of the paper, we consider a function / in the form 

/(e, 6) = exp (A | ^ | /3 sgn(e)) . (ft a l e>0 + l e =o + «&l e <o) , (2.5) 

where /3 > is the price impact exponent, A > is the temporary price impact factor, 
< 1, and K a > 1 are the bid and ask spread parameters. The impact of liquidity modelled 



in (2.4) is like a transaction cost combining nonlinearity and proportionality effects. The 
nonlinear costs come from the dependence of the function / on e, but also on 9. On the 
other hand, this transaction cost function / can be determined implicitly from the impact 
of a market order placed by a large trader in a limit order book, as explained in [15] , [T9] 



or |18| . Moreover, the dependance of / in 9 in (2.5) means that rapid trading has a larger 



temporary price impact than slower trading. Such kind of assumption is also made in the 



seminal paper [2], and reflects stylized facts on limit order books. The form (2.5) was 
suggested in several empirical studies, see [12], [T7], [3], and used also in [8], [TT] . 

Cash holdings. We assume a zero risk-free return, so that the cash holdings are constant 
between two trading times: 

X t = X Tn , r n <t<T n +i, n>0. (2.6) 

When a discrete trading Al^ = Cn+i occurs at time t = r n +i, this results in a variation of 
the cash amount given by AX t := X t — X t - = —AY t .Q(AY t , P tl @ t -) due to the illiquidity 
effects. In other words, we have 

X r n+1 = X r -^ - Cn+lPr n+1 f(Cn+l,r n+1 - T n ), U > 0. (2.7) 

Remark 2.1 Notice that since /(e, 0) = if e < and /(e, 0) = oc if e > 0, an immediate 
sale does not increase the cash holdings, i.e. X Tn+1 = X - = X Tn1 while an immediate 

n+l 

purchase leads to a bankruptcy i.e. X Tn+1 — — oc. 



Liquidation value and solvency constraint. The solvency constraint is a key issue in 
portfolio choice problem. The point is to define in an economically meaningful way what 
is the portfolio value of a position in cash and stocks. In our context, we first impose a 
no-short selling constraint on the trading strategies, i.e. 

Y t > 0, < t < T. 

Next, we introduce the liquidation function L(x,y,p,9) representing the value that an 
investor would obtain by liquidating immediately his stock position y by a single block 
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trade, when the pre-trade price is p and the time lag from the last order is 9. It is defined 
onRxR+x (0, oo) x [0, T] by 

L(x,y,p, 9) = x + ypf(-y,0), 

and we constrain the portfolio's liquidative value to satisfy the solvency criterion: 

L{X u Y u P u ® t ) > 0, 0<£<T. 

We then naturally introduce the solvency region: 

S = {(z, 9) = (x, y,p,9) ERx R + x (0, oo) x [0, T] :L(z,9)>0}. 

and we denote its boundary and its closure by 

dS = d y S U d L S and 5 = 5U dS. 

where 

d y S = {(z,0) = (x,y,p,0) ERx R + x (0, oo) x [0, T] : y = and x = L(z,0) > 0} , 
d L S = {(z, 9) = (x, y,p,9) ERx R + x (0, oo) x [0, T] : L(z,9) = 0}. 

In the sequel, we also introduce the corner lines in dS : 

D = {(0,0)} x (0,oc) x [0,T] = d y Snd L S. 

Admissible trading strategies. Given (£, z, 9) E [0, T] x <S, we say that the impulse control 
strategy a = (r n , ( n )n>o is admissible, denoted by a E A(t, z, 9), if to = t — 9, r n > £, n > 1, 
and the process {(^ S ,6 S ) = (X s F s , P s , 9 S ), t < s < T} solution to d27Tb - d^T^ - d27^ - dzTsb - 



( |2.7[ ), with an initial state (Z t -,@ t -) = (z,9) (and the convention that (Zt,@t) — ( z i@) if 
ri > £), satisfies (Z s , S ) E [0,T] x 5 for all s E [£, T]. As usual, to alleviate notations, we 
omit the dependence of (Z, O) in (£, z, 9, a), when there is no ambiguity. 

Portfolio liquidation problem. We consider a utility function U from R + into R, strictly 
increasing, concave and w.l.o.g. U(0) = 0, and s.t. there exists K > 0, 7 E [0, 1) : 

J7(w) < Kid 7 , Vw > 0. 

The problem of optimal portfolio liquidation is formulated as 

v(t,z,9) = 



sup E [U(X T )] , (t, z, 6) € [0, T] x 5, 



(2i 



where Ai(t,z,9) — {a E A(t,z,0) : It = 0}. As observed in [IT], one can shift the 
terminal liquidation constraint in Ai(t, z, 9) to a terminal liquidation utility by considering 
the function Ul defined on S by: 

U L (z,9) = U(L(z,9)), (z,9)€S. 

Then, problem (|2.8|) is written equivalently in 



v(t,z,9) = sup E 

aeA(t,z,0) 



u L (z T ,e T ) , (t,z,0) e [0,T] xS. 



(2.9) 
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2.2 PDE characterization 



The dynamic programming Hamilton-Jacobi-Bellman (HJB) equation corresponding to the 
stochastic control problem (2.8) is a quasi- variational inequality written as 

dv — 
min [ - — - Cv , v - Hv] = 0, on [0, T) x S, (2.10) 

together with the relaxed terminal condition 

mm[v-U L ,v-Hv] = 0, on {T} x S. (2.11) 

Here, C is the infinitesimal generator associated to the process (Z = (X, Y, P),0) in a 
no-trading period: 

dip n dip 1 n oflV 

% is the impulse operator defined by 

H<p(t, z, 9) = sup <p(t, T(z, 9, e), 0), (i, z, 0) E [0, T] x 5, 

eeC(z,0) 

r is the impulse transaction function defined from 5xR into R x R x (0, oo): 

T(z,e,e) = (x-epf(e,0),y + e,p), z = (x, y,p) e 5, e G M, 

and C{z 1 9) the set of admissible transactions : 

C(*,0) = {eGR : (r(z,0,e),o) E 5}. 

By standard arguments, we derive the constrained viscosity solution property of the 
value function v to (2.10)-(2.11 ). However, in order to have a complete characterization of 
the value function via its HJB equation, we need a uniqueness result. Unfortunately, in our 
model, it seems not possible to get such result, at least by classical arguments since there 
is no strict supersolution to (2.10). In [11], the authors prove a weaker characterization of 
the value function in terms of minimal solution to its HJB equation. They also consider 
a small variation of the original model by adding a fixed transaction fee e > at each 
trading. This means that given a trading strategy a = (r n ,(^ n ) n >o, the controlled state 
process (Z — (X, Y, P), 0) jumps now at time r n+ i, by: 



(z Tn+1 ,e Tn+1 ) = [v e {z e T - Cn+i),o , (2.12) 

where T £ is the function defined onKxM + x (0, oo) x [0, T]xR into iU {-00} xKx(0, 00) 
by: 

T e (z,0,e) = r(z,6,e)-(e,0,0) = (x - epf(e,0) - e, y + e,p) , 

for z — (x,y,p). The dynamics of (Z,Q) between trading dates is given as before. We 
introduce a modified liquidation function L £ defined by: 

L £ (z,6) = max[x,L(z,0)-e], (z,0) = (x,y,p,9) € R x R + x (0,oo) x [0,T]. 
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The interpretation of this modified liquidation function is the following. Due to the presence 
of the transaction fee at each trading, it may be advantageous for the investor not to 
liquidate his position in stock shares (which would give him L(z,9) — e), and rather bin his 
stock shares, by keeping only his cash amount (which would give him x). Hence, the investor 
chooses the best of these two possibilities, which induces a liquidation value L £ (z,9). 

The corresponding solvency region S £ C S with its closure S £ = S £ U dS £l and boundary 
dS £ = d y S £ U 3lS £ are given by: 

S £ = | (z, 9) = (x, y,p,9) E R x R + x (0, oo) x [0, T] : y>0 and L e (z, 9) > o}, 
d y S £ = | (z, 9) = (x, y,p,9) Elx R + x (0, oo) x [0, T] : y = and L £ (z,6>)>o}, 
d L S £ = | (z, 9) = (x, y,p,9) Elx R + x (0, oo) x R + : L £ (z,6>) = o}. 

The set of admissible trading strategies is defined as follows: given (£, z, 9) E [0, T] x5 e , 
we say that the impulse control a is admissible, denoted by a E *A e (£, 0), if To = t — 0, r n 
> t, n > 1, and the controlled state process (Z £ , 9) solution to (|^-(|^-(|^-(^ 
with an initial state (Z^_,Q t -) = (z,9) (and the convention that (Z £ ,@t) = (z,9) if t\ > 
£), satisfies (Zf,0 s ) E [0,T] x <S £ for all s E [t, T]. Here, we stress the dependence of Z £ = 
(X £ , Y, P) in e appearing in the transaction function T £ , and we notice that it affects only 
the cash component. 

The liquidation utility function in this model with fixed transaction fee e is defined on 
S £ by Ul £ (z,9) = U(L £ (z,9)), and the associated optimal portfolio liquidation problem is 
defined via its value function by: 

v e (t,z,0) = sup E[U Le (Z^e T )], (t,z,6)€[0,T]xS e . (2.13) 

aeA £ (t,z,0) 



The dynamic programming equation associated to the control problem (2.13) is 

dv, 



mm 



on [0,T)x5 £ , (2.14) 



- ~ £>V £ , V £ ~ U £ V £ 

mm[v £ -U Le ,v £ -n £ v £ ] = on {T}x5 e , (2.15) 



where % £ is the impulse operator defined by 

H e w(t,z,6) = sup w(t,r e (*,0,e),O), (£, z, 9) E [0, T] x <S £ , 

for any locally bounded function w on [0,T] x <S £ , with the convention that % e w(t,z,Q) 
— — oo when C £ (z,9) = 0, and the set of admissible transactions in the model with fixed 
transaction fee is: 

C e (z,0) = {eeR:(r £ ( 2 ,^),0)e5 £ }, (z, 6) G S £ . 

We recall from [IT] that v £ is in the set £/([0, T] x S e ) of functions satisfying the growth 
condition: 

g([0,T]xS £ ) = L:[0,T]x5 ( ^Ms.t. sup < oo}. 



[0 
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Remark 2.2 The function z v £ (t,z,0) is strictly increasing in the argument of cash 
holdings x, for (z = (x,y,p),0) E <S £ , and fixed t E [0, T]. Indeed, for x < x f , and z 
= (x,y,p), z f — (x',y,p), any strategy a E A £ (t,z,9) with corresponding state process 
(Z s = (X s , y s , P s ), O s ) s >t, is also in A £ (t, z f ,9), and leads to an associated state process 
(Zg = (X s + x r — x 1 Y Sl P s ) 1 Q s ) s >f Using the fact that the utility function is strictly 
increasing, we deduce that v £ (t,x,y,p,0) < v £ (t, x f , y,p, 0). Moreover, the function z —> 
v £ (t, z, 0) is nondecreasing in the argument of number of shares y. Indeed, fix z = (x, y,£>), 
and z' — (x,y f ,p) with y < y' . Given any arbitrary a = (r n Xn)n £ A £ (t,z,0), consider 
the strategy ol — (t^,^), starting from (x,y',p) at time t, which consists in trading again 
immediately at time t by selling y' — y shares (which does not change the cash holdings, see 
Remark |2.1[ ), and then follow the same strategy than a. The corresponding state process 
satisfies {Z' s1 & s ) — (Z s , @ s ) a.s. for s > t, and in particular ol E A £ (t, z\ 0), together with 
E[[/l £ (Z^, O^)] = K[Ul £ (Zt, St)] < v(t, z f , 9). Since a is arbitrary in A £ (t, z, 0), this shows 
that v(t,x,y,p,0) < v(t,x,y',p,0). 

In the sequel, we shall denote by <5+([0, T] x S £ ) the set of functions tp in <5([0, T] x S £ ) 
such that (p(t, x, y,p, 0) is strictly increasing in x and nondecreasing in y. 

Remark 2.3 Fix t E [0, T]. For = 0, and z = (x, y,p) s.t. (z, 0) E <S £ , the set of admissible 
transactions C e (z, 0) = [— y, 0] (and r e (z, 0, e) = (x — 5, y + e,p) for e E C £ (z, 0)) if x > e, 
and is empty otherwise. Thus, % £ w{t, z, 0) = sup eG [_ y ] ^(^ x — 6, y + e,p, 0) if x > £, and 
is equal to — oc otherwise. This implies in particular that 

H £ w(t,z,0) < w(t,z,0), (2.16) 



for any w E £7+([0, T] x <S £ ), which is the case of v £ (see Remark 2.2). Therefore, due to the 



market impact function / in (2.5) penalizing rapid trades, it is not optimal to trade again 



immediately right after some trade, i.e. the optimal trading times are strictly increasing. 

A main result in [llj is to provide a unique PDE characterization of the value functions 
v £ , e > 0, and to prove that the sequence (v £ ) £ converges to the original value function v 
as e goes to zero. 

Theorem 2.1 (1) The sequence (v £ ) £ is nonincreasing, and converges pointwise on [0,T] x 
(S \ OlS) towards v as e goes to zero, with v £ < v. 

(2) For any e > 0, the value function v £ is continuous on [0,T) x S £ , and is the unique (in 
[0,T) x S £ ) constrained viscosity solution to (2.14)-( 2.15| ), satisfying the growth condition 
in £/([0, T] x S £ ), and the boundary condition: 

lim v £ (t f ,z',9') = v e (t,z,0) 
(tf,*',0')-KM,0) 

= C7(0), V(t,^(0,0,ri^)G[0,T]xD o . (2.17) 
The rest of this paper is devoted to the numerical analysis and resolution of the QVI 



(2.14)-(2.15), which then provides an ^-approximation of the original optimal portfolio 



liquidation problem (2.8). 
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3 Time discretization and convergence analysis 



In this section, we fix s > 0, and we study time discretization of the QVI (2.14)-(2.15) 
characterizing the value function v £ . For a time discretization step h > on the interval 
[0,T], let us consider the following approximation scheme: 

S h (t,z,9,v h (t,z,6),v h ) 



0, (t,z,0)€[0,T\ 
where S h : [0, T] x S £ x R x (? + ([0, T] x S e ) -»• R is defined by 
S h (t,z,0,r,<p) 



(3.1) 
(3.2) 




r-E[ V (t + fc ) z!j5f ) eW j 



H e <p(t,z,6) 



r-E[ip(T,Z^ z ,G°/> d )] , r-U e i P {t,z,e) 
r-U Le (z,9) , r-H e ip(t,z,8) 



if t e [0, r - /i] 
if t € (T - /i, T) 
if t = T. 



Here, (Z°'*> z , 0°'*^) denotes the state process starting from (z, 0) at time t, and without 
any impulse control strategy: it is given by 



( Z 0,M @ 0, t)( 



with the solution to (2.3) starting from p at time t. Notice that (3.1) is formulated as 
a backward scheme for the solution v h through: 

v h (T,z,9) = max[U Lc (z,0) , H e v h (T, z,9)], 



v h (t,z,9) 



max 



E[v h (t + h, Z%£, 9 + h)} , U e v h (t, z, 9) 



(3.3) 

0<t<T-h, (3.4) 



and v h (t,z,9) — v h (T — h, z, 9) for T — h < t < T. This approximation scheme seems a 
priori implicit due to the nonlocal obstacle term % e . This is typically the case in impulse 
control problems, and the usual way (see e.g. [§], [H]) to circumvent this problem is to 
iterate the scheme by considering a sequence of optimal stopping problems: 

v h ' n+1 (T,z,0) 



max[U Le (z,e) , H e v h ' n {T,z,e)], 



v h > n+1 (t,z,6) 



max 



E[v h ' n+1 (t + h, Z^' h z , 9t + h)] , Hev h > n (t, z, 9) 



starting from v h >° = E[U Ls (Z^' z , Of*' )]. Here, we shall make the numerical scheme ( |3.l[ ) 
explicit, i.e. without iteration, by taking effect of the state variable 9 in our model. Recall 
indeed from Remark 2.3 that it is not optimal to trade again immediately right after some 
trade. Thus, for v h E £?+([0,T] x <S £ ), and any (V,0) E <S £ , we have from (2.16) and 

v h (T,z\0) = U Le (z\0) 
v h (t,z?,0) = E[v\t + h,$tf,h)]. 



Therefore, by using again the definition of T~L £ in the relations ( |3.3| )-( [3^4| ), we see that the 
scheme (|3.1|) is written equivalently as an explicit backward scheme: 

(3.5) 

E [v h (t + h, Z^' z , 9 + h)], sup E [v h (t + h, Z° t '!f e ,h)]], (3.6) 

eec e (z,e) J 



v h (T,z,9) = max[U Le (z,9) , U e U Le (z,e)], 
v h (t,z,9) 



max 
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for < t < T - h, and v h (t, z, 9) = v h (T - h, z,0) for T - h < t < T, where we denote z% 



T £ {z,Q,e) in (3.6) to alleviate notations. Notice that at this stage, this approximation 



scheme is not yet fully implementable since it requires an approximation method for the 



expectations arising in (3.6). This is the concern of the next section. 



We focus now on the convergence (when h goes to zero) of the solution v h to (3.1) 
towards the value function v £ solution to ( |2.14 )-( 2.15 ). Following [3], we have to show 
that the scheme S h in (3.2) satisfies monotonicity, stability and consistency properties. As 



usual, the monotonicity property follows directly from the definition (3.2) of the scheme 



Proposition 3.1 (Monotonicity) 
For all h > 0, (t, z, 9) G [0, T] x S e , r G 



and (p, ifj G £?+([0,T] x S e ) s.t. <p < ifj, we have 
S h (t,z,9,r,cp) > S h (t,z,e,r^). 



We next prove the stability property. 

Proposition 3.2 (Stability) 
For all h > 0, there exists a unique solution v h G 

(v h )h is uniformly bounded in Q([0,T] x S £ ): there exists w G C/([0 5 T] x S £ ) s.t. 
for all h > 0. 

Proof. The uniqueness of a solution G 



T] x S £ ) to (3.1) ; and the sequence 



< \w\ 



, T] x S £ ) to (3.1) follows from the explicit 



backward scheme (3.5)-(3.6). For t G [0,T], denote by N^h l ne integer part of (T — t)/h, 
and T tj h = {tk — t + fc/i, k = 0, . . . , N t ^} the partition of the interval [£, T] with time step 
h. For (t,z,9) G [0,T] x <S £ , we denote by A £ (t,z,9) the subset of elements a = (r n ,C n ) n 
in A £ (t,z,9) such that the trading times r n are valued in T^. Let us then consider the 
impulse control problem 



v\t,z,0) 



sup E[U Le (Z^,e T ) 

aeA*(t,z,0) 



(t,z,9) G [0,T] xS s 



(3.7) 



It is clear from the representation (3.7) that for all h > 0, < v h < v £l which shows that 
the sequence (v h )h is uniformly bounded in £/([0, T] x S £ ). Moreover, similarly as for v £l and 
by the same arguments as in Remark 2.2, we see that v h (t,z,0) is strictly increasing in x 



and nondecreasing in y for (z, 0) = (x, 0) G S £ . Finally, we observe that the numerical 
scheme (3.1) is the dynamic programming equation satisfied by the value function v h . This 



proves the required stability result. 

We now move on the consistency property. 

Proposition 3.3 (Consistency) 

(i) For all (t,z,9) G [0,T) x S e and cf) G C^ 2 ([0,T) x S € ), we have 



□ 



lim sup min < 

(h,t' ,z ,0')^-(O,t,z,0) 

(V,z / ,6' / )G[0,T)x<S e 



■E 



+ h > Z t> +h > @ t' +h > 



< min 



-{( 



dt 



h 



Ccf>)(t,z,6), 



; <i>)(t,z,e') 



!>)(t,z,0)} 



(3.J 



li 



and 



lim inf 



mm < 



(h,t' ,z ,0')^-(O,t,z,0) 

It' ,z ,0')e[o,T)xS e 



h 



I I J X 



> 



mm 



{ ( - ^ - £4) (t, z,e), (0 - u e <j>) (t, z,o)} 

(ii) For all (z,6) € <S e and <j> € C 1>2 ([0,r] x 5 e ), we faroe 



(3.9) 



and 



lim sup 

(V ,z ,6' ')^(T,z,8) 

(t f ,z ,0')e[0,T)xS € 


min , z , 9 ) 


-U Le (z,0'), 




< 


min |^(T, z, 0) - 


U Le (z,9),^ 


-U e $)(T,z,9)} 


lim inf 

(V ,z ,6' ')^(T,z,8) 

(t' ,z' ,e')e[o,T)xS e 


min , z , ) 


-U Le (z,9 f ), 


(<t>-U^){t,z,9') 


> 


min | (V(T, z, 0) 




(cj>-U e ^){T,z,9)} 



(3.10) 



(3.11) 

Proof. The arguments are standard, and can be adapted e.g. from [6] or [7J. We sketch 
the proof, and only show the inequality (3.8) since the other ones are derived similarly. Fix 
t E [0, T). Since the minimum of two upper-semicontinous (use) functions is also use and 
using the caracterization of use functions, we have 



t\z\e') 



E 



lim sup min { ( — - 

z ,6 f )e[0,T)xS e 

< lim sup min < lim sup ((f) — T-L e (f)j (t ,z ,9 ), 



, 7 ry0,t ,Z r^0,t 



(h,t ,z ,0 )-KO,7j,z,0) 

(tV,0')e[o,T)x<s e 



(t'V,0")e[O,T)xS e 



// // ^// 



z ,9 ) — E 



lim sup 

'V)-Ko,i 

",e")e[o,i 
< min < lim sup 

^ /, ./ / „/x 



(m'VV')^(o,*V/) 

(t",z",e")e[o,T)xs e 



(h,t ,z ,0 )->(0,t,z,9) 

(t' ,z' ,9')<E[0,T)xSe 



I I J X 



HM(t,z,0), 



#', Z >')-E[#' + M^,e^)] 

limsup b > 

h ) 



(h,t' ,z' ,e f )^(o,t,z,o) 
(tV,0')E[O,T)x«S e 



< min{0(t,z,0) -H € <l>(t,z,0) 



lim sup 

(^,t / ,z / ,6» / )^(0,t,2,6») 

(tV,0')E[o,T)x<s e 



0(* ,2 , 



E 



(3.12) 
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where the last inequality follows from the continuity of (j) and the lower semicontinuity of 

a i! z 1 i! r 

% £ . Moreover, by Ito's formula applied to 0(s, Z s ' ' , S ' ' ), and standard arguments of 

localization to remove in expectation the stochastic integral, we get 



lim sup 

(h/ ,z' ,0')->(0,t,z,6) 

it' ,z' ,6')(E[Q,T)xS e 



(f)(t ,Z ,1 



E 



■ 7 ryO.t ,Z r^O.t 



t'+h ' 



h 



+ 



Substituting into (3.12), we obtain the desired inequality (3.8). 



□ 



Since the numerical scheme (3.1) is monotone, stable and consistent, we can follow the 
viscosity solutions arguments as in [4J to prove the convergence of v h to v e , by relying on 
the PDE characterization of v £ in Theorem |2.1| (2), and the strong comparison principle for 
(|27l4|)-(|27l5|) proven in [IT]. 



Theorem 3.1 (Convergence) The solution v h of the numerical scheme (3.1) converges 
locally uniformly to v £ on [0, T) x S e . 



Proof. Let v e and v e be defined on [0,T] x S £ by 

vl(t,z,0) — lim sup 



(h,t' ,z' ,e' ')^(o,t, z ,e] 
(t f ,z' ,e')£[0,T)xS t 



Ve(t,Z,0) 



lim inf 



(h,t' ,z' ,e' ')^(o,t, z ,e) 

(t f ,z ,e')e[o,T)xS £ 



v h (t\z\e') 

v h {t,z^) 



We first see that % and are respectively viscosity subsolution and supersolution of (2.14)- 
( |2. 15 ) . These viscosity properties follow indeed, by standard arguments as in [3] (see also 
[6] or [7] for impulse control problems), from the monotonicity, stability and consistency 
properties. Details can be obtained upon request to the authors. Moreover, from (3.7), we 
have the inequality: £7(0) < v h < v £ , which implies by (2.17): 

liminf vj(t \z\e') = £7(0) = l%(t,z,0), V (t, z, 9) E [0, T] x Dq (3.13) 

(t' \z ,e' ')^(t,z,e) 
(t' ,z ,o')e[Q,T)xS e 

Thus, by using the strong comparison principle for (2.14)-( 2.15[ ) stated in Theorem 5.2 [11], 
we deduce that W e < on [0, T] x S e and so v~ e = = v e on [0, T] x S € . This proves the 
required convergence result. □ 



4 Numerical Algorithm 

Let us consider a time step h = T/m, m E N \ {0}, and denote by T m = {t{ = ih,i — 
0, . . . , m} the regular grid over the interval [0, T]. We recall from the previous section that 
the time discretization of step h for the QVI ( 2. 14| )-( 2. 15 ) leads to the convergent explicit 
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backward scheme: 



v h (t m ,z,0) 



max 



U Le (z,0) , 

sup v h (t m ,F £ (z,6,e),0) 

eeC e (z,0) 



v h (t h z,9) 



E[v h (t 
max 



l+1 ,z^f,e + h)] 



E[v h (t i+1 ,Z^; z ,e + h)] , 

sup i; ft (*i,r e (z,0,e),O)] 
eec £ (z,e) 



if e = o 
if e > o, 

if = 
if > 



(4.1) 



(4.2) 



for i = 0, . . . , m — 1, (2 = (x,y,p),0) E <S £ . Recall that the variable represents the time 
lag between the current time £ and the last trade. Thus, it suffices to consider at each time 
step U of T m , a discretization for 9 valued in the time grid 

Ti, m {6j=jh, j = 0, i = 0,...,ra. 

On the other hand, the above scheme involves nonlocal terms in the variable 2 for the 
solution ^ in relation with the supremum over e E C £ (z 1 9) and the expectations in (4.1)- 
(4.2), and thus the practical implementation requires a discretization for the state variable 
z, together with an interpolation. For any 0j E T^ m , let us denote by 

Z j = {z = (x,y,p) E R x R + x R + : (z,^-) E <S £ }- 

For the discretization of the state variable z E Z\ and since <S £ is unbounded, we first 
localize the domain by setting Z\ oc = Z J n ([%m,%aj x [VmimVmax] x bmm^mox]), 



< ^ 



m 



< u 



oc 

mm ^ Vmaxi 



< Pmin < Pmax are fixed constants, and then 



where x min 
define the regular grid: 

K = {z = (x,y,p) eX n xY n xF n :(z,9j) eS £ }. 

where X n is the uniform grid on [x m i ni x max ] of step Xmax ~ x ™ in ; and similarly for Y n , 



Optimal quantization method. Let us now describe the numerical procedure for com- 
puting the expectations arising in d42b. Recalling that Z '*' 2 = (x,y,P*P), this involves 
only the expectation with respect to the price process, assumed here to follow a Black- 
Scholes model (2.3). We shall then use an optimal quantization for the standard normal 
random variable J7, which consists in approximating the distribution of U by the discrete 
law of a random variable U of support (uk)i<k<N £ ^ N \ and defined as the projection 
of U on the grid (uk)i<k<N according to the closest neighbour. The grid (uk)i<k<N is 
optimized in order to minimize the distorsion error, i.e. the quadratic norm between U 
and U. This optimal grid and the associated weights (7Tk)i<k<N are downloaded from the 
website: http://87.106.220.249/n01. We refer to the survey article [16] for more details on 
the theoretical and computational aspects of optimal quantization methods. From ( |4.2[ ), 
we have to compute at any time step U E T m , and for any 9j E T^ m , z = (x,y,p) E Z^, 
expectations in the form: 

£ h {t i ,z,e j ) := E[v h (U + h,Z%*g,0 j + h)] 

2 

= E[v h (ti + h,x,y,pexp((b- %-)h + aVhU),0j + h)], 
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that we approximate by 

1 N a 2 
£ h (ti,z,9j) ~ — J2^kV h (t i + h,x,y,pexp((b-^)h + aVhu k ),9 j + h). (4.3) 

k=i 



Interpolation procedure. Notice that for implementing recursively in (4.2) this quan- 
tization method, we need to compute z —> v h (ti, z,9j) on Z J given the known values of 
v h (ti, z,9j) on 7? n . This approximation is achieved as follows: 

• if z G 2j oc1 then we use a linear interpolation of v h {ti, z,0j) with respect to the closest 
neighbours in l? n . 

• if z = (x, y,p) ^ Zi oc , we use the growth condition satisfied by the value function: 

v h (t i ,z,9 i ) ~ v h (t l ,z,9 1 ) [ X + yP } 1 

where z = (x,y,p) is the projection of z (according to the closest neighbour) on the 
grid Z 3 n . 

Algorithm description. In summary, our numerical scheme provides an algorithm for 
computing approximations v h of the value function, and ( h of the optimal trading strategy 
at each time step U G T m , and each point (z, 9) of the grid (Z n x T^ m ) D S £ . The parameters 
in the algorithm are: 

- T the maturity 

- b and a the Black and Scholes parameters of the stock price 



- A the impact parameter, f3 the impact exponent in the market impact function (2.5) 

- K a , Kb the spread parameters in percent, e the transactions costs fee 

- We take by default a CRRA utility function: U{x) = x 1 

G R, < y min < ymax, < p min < p ma x, the boundaries of the localized 

domain 

- m number of steps in time discretization, n the number of steps in space discretization 

- N number of points for optimal quantization of the normal law, M number of points used 
in the static supremum in e 

The algorithm is described explicitly in backward induction as follows: 
► Initialization step at time t m — T: 

• (s:0) For j = 0, set v h (t m , z, 0) = Ul £ (z,0), (> h {t m ,z,{)) = on Z^, and interpolate 
v h (t m ,z,0) on Z°. 

• (s:j) For j = 1, . . . ,ra, 

— for z G Z^, compute u := sup L r L £ (r £ (z, 9j, e), 0) and denote by e the argu- 

e6C e (*,0j) 

ment maximum: 

— if v > UL e ( z i@j)i then set v h (t mi z,0j) — v and ( h (t mi z,0j) — e, 
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— else set v h (t mi z,Oj) = UL £ {z,Oj), and ( h (t mi z,0j) — 0. 

- Interpolate z — )► v h (t mi z,9j) on ZK 

► From time step t{+i to U, i = m — 1, . . . , 0: 



fs.-Oj For j = 0, compute £ h (ti,z,0) from (4.3) and (s:l) of time step t^+i, and set 



^(ti,^,0) = *,0), ^(tij^O) = on interpolate ^(^,z,0) on Z°. 

(s:j) For j = 1, . . . ,i, 



— for 2 E Z^, compute £ h (ti,z,9j) from (4.3) and (s:j+l) of time step i; : 



sup ^(ti, r e (z, 9j, e), 0) from fs/Oj, and denote by e the argument maximum: 

eeC e (z,0j) 

— if v > £ h (ti, z,0j), then set v h {ti,z,0j) — v, and (> h {ti,z,0j) — e, 

— else set v h (ti,z,9j) = £ h (ti,z,9j), and ( h (ti,z,9j) = 0. 

— Interpolate z v h {ti,z,0j) on Z J . 

5 Numerical Results 
5.1 Procedure 

For each of the numerical tests, we used the same procedure consisting in the following 
steps: 

(1) Set the parameters according to the parameter table described in the first subsection 
of each test 

(2) Compute and save the grids representing value function and optimal policy according 
to the optimal liquidation algorithm 

(3) Generate Q paths for the stock price process following a geometrical Brownian motion: 
we choose parameters b and a that allows us to observe several empirical facts on the 
performance and the behavior of optimal liquidation strategy. These parameters can also 
be estimated from historical observations on real data by standard statistical methods. 

(4) Consider the portfolio made of Xq dollars and Yq shares of risky asset 

(5) For each price path realization, update the portfolio along time and price path accord- 
ingly to the policy computed in the second step 

(6) Save each optimal liquidation realization 

(7) Compute statistics 

In the sequel, we shall use the following quantities as descriptive statistics: 
• The performance of the i-th realization of the optimal strategy is defined by 

opt X + Y P 
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where {Z^ h , 0y- ) is the state process, starting at date at (Xo, Yo, Po, 0), evolv- 
ing under the z-th price realization and the optimal control a opt . This quantity can be 
interpreted as the ratio between the cash obtained from the optimal liquidation strategy 

(i) 

and the ideal Merton liquidation. We define in the same way the quantities ^ naive and 
L^l-r _ respectively associated with the controls a naive and a uni ^ orm of the naive and 

uniform 1 J 

uniform strategy, refereed as benchmark strategies. Recall that the naive strategy consists 
in liquidating the whole portfolio in one block at the last date, and the uniform strategy 
consists in liquidating the same quantity of asset at each predefined date until the last date. 
Notice that the score 1 corresponds to the strategy, which consists in liquidating the whole 
portfolio immediately in an ideal Merton market. 

When denoting by Q the number of paths of our simulation, we define: 



• The mean utility V Tj^2 u ( L ^ ] 



1 Q 

• The mean performance L=— 



• The standard deviation of the strategy <r 



^S (i,,,)2 - £i 



Here the dot . stands for opt, naive or uniform. We will also compute the third and fourth 
standardized moments for the series 



5.2 Test 1: A toy example 

The goal of this test is to show the main characteristics of our results. We choose a set 
of parameters that is unrealistic but that has the advantage of emphasizing the typical 
behavior of the optimal liquidation strategy. 

Parameters We choose the set of parameters shown in table [IJ 

Execution statistics The results were computed using Intel® Core 2 Duo at 2.93Ghz 
CPU with 2.98 Go of RAM. Statistics are shown in tabled 

Shape of policy In this paragraph we plotted the shape of the policy sliced in the plane 
(x,y), i.e. the (cash, shares) plane, for a fixed (t,0,p) (figure [I]). The color of the map 
at (#05 2/o ) on the graph represents the action one has to take when reaching the state 
(t, 0, #o, 2/0 j We can see three zones: a buy zone (denoted BUY on the graph), a sell 
zone (denoted SELL on the graph) and a no trade zone (denoted NT on the graph). Note 
that the bottom left zone on the graph is outside the domain S. These results have the 
intuitive financial interpretation: when x is big and y is small, the investor has enough cash 
to buy shares of the risky asset and tries to profit from an increased exposure. When y is 
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Parameter 


Value 


Parameter 


Value 


Maturity 


1 year 


X 


2000 


A 


5.00E-07 


Y 


2500 





0.5 


Po 


5.0 


7 


0.5 


%min 


-30000 


HA 


1.01 


%max 


80000 


K B 


0.99 


Umin 





e 


0.001 


Umax 


5000 


b 


0.1 


Pmin 





a 


0.5 


Pmax 


20 






rn 


40 






n 


20 






N 


100 






Q 


10 5 



Table 1: Test 1: parameters 



Quantity 


Evaluation 


Time Elapsed for grid computation in seconds 


7520 


Number Of Available Processors 


2 


Estimated Memory Used (Upper bound) 


953MB 


Time Elapsed for statistics Computation in seconds 


21 



Table 2: Test 1: Execution statistics 



large and x is small, the investor has to reduce exposure to match the terminal liquidation 
constraint. 

We also plotted the shape of the policy sliced in the plane i.e. the (shares,price) 

plane, for a fixed (t, 0, x) (figure [2]). As before, the color of the map at (yo,po) on the graph 
represents the action one has to take when reaching the state (£, 0, x, yo,po). Again, we can 
distinguish the three zones: buy, sell and no trade. 

Shape of value function Figure [5] shows the value function sliced in the (x, y) plane. 
This figure is a typical pattern of the value function. Recall from Proposition 3.1 in 
the following Merton theoretical bound for the value function: 

v(t,z,0) < v M (t,x,y,p) = e p(T_t) (x + ^) 7 , with p = 7^—^- 

1 — 7 la* 

In the figure [4] we plotted the difference between the value function and this theoretical 
bound. We observe that this difference is increasing with the number of shares, and de- 
creasing with the cash. This result is interpreted as follows: the price impact increases with 
the number of shares, but this can be reduced by the liquidation strategy whose efficiency 
is greater if the investor can sustain bigger cash variations. 
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Policy sliced in the (cash,shares) plane 




5 # c# & <P ^ -s*> 4^ ^ ^ ^ 

J0 5 # 4> N tS> rfi ng> $ tfP <^ «^ ^ A 



Figure 1: Test 1: Typical shape of the policy sliced in XY near date 



Policy sliced in the (price,shares) plane 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 

PRICE 



Figure 2: Test 1: Typical shape of the policy sliced in YP 
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Value Function typical shape sliced in XY 




Figure 3: Test 1: Typical shape of the value function sliced in XY 



Difference between value function and Merton liquidation value 




Figure 4: Test 1: Difference between value function and Merton theoretical bound 
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5.3 Test 2: Short term liquidation 

The goal of this test is to show the behavior of the algorithm on a realistic set of parameters 
and real data. We used Reuters™ data fed by OneTick™ TimeSeries Database. We used 
the spot prices (Best Bid and Best Ask) for the week starting 04/19/2010 on BNP.PA. We 
computed mid-price that is the middle between best bid and best ask price. We choose the 
impact parameter A in order to penalize by approximately 1% the immediate liquidation 
of the whole portfolio compared to Merton liquidation. In other words, we take A so that: 

A&-0.01. 
i T 

Parameters We computed the strategy with parameters shown in table [3j 



Parameter 


Value 


Parameter 


Value 


Maturity 


1 Day 


X 


20000 


A 


5.00E-04 


Yq 


2500 


P 


0.2 


Po 


52.0 


7 


0.5 




-30000 


K A 


1.0001 


%max 


200000 




0.9999 


Umin 





e 


0.001 


Umax 


5000 


b 


0.005 


Pmin 


50.0 


a 


0.25 


Pmax 


54.0 






rn 


30 






n 


40 






N 


100 






Q 


10 5 



Table 3: Test 2: Parameters 



Execution statistics We obtained the results using; Intel® Core 2 Duo at 2.93Ghz CPU 
with 2.98 Go of RAM, the computations statistics are gathered in table [1J 



Quantity Evaluation 

Time Elapsed for grid computation in seconds 8123 

Number Of Available Processors 2 

Estimated Memory Used (Upper bound) 573MB 



Table 4: Test 2: Execution statistics 

Performance Analysis We computed the mean utility and the first four moments of 
the optimal strategy and the two benchmark strategies in table [5] and plotted the empirical 
distribution of performance in figure [5j It is remarkable that the optimal strategy gives an 
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empirical performance that is above the immediate liquidation at date in the Merton ideal 
market. This is due to the fact that the optimal strategy has an opportunistic behavior, 
as the decisions are based on the price level, and so profit from the 'detection" of some 
favorable price conditions. Indeed, an optimal trading strategy is embedded with the 
optimal liquidation: in this example, this feature not only compensates the trading costs, 
but also provides an extra performance compared to an ideal immediate liquidation at date 
0. Still, the Merton case is a theoretical upper bound in the following sense: the optimal 
value function with trading costs is below the optimal value function without trading costs, 
recall the figure [4j As expected, the empirical distribution is between the distributions of 
the two other benchmark strategies. We also notice that the optimal strategy outperforms 
the two others by approximatively 0.25% in utility and in performance. 



Strategy 


Utility V 


Mean L 


Standard Dev. 


Skewness 


Kurtosis 


Naive 


0.99993 


0.99986 


0.00429 


0.94584 


4.68592 


Uniform 


0.99994 


0.99988 


0.00240 


0.42788 


3.34397 


Optimal 


1.00116 


1.00233 


0.00436 


1.03892 


4.89161 



Table 5: Test 2: Utility and first four moments for the optimal strategy and the two benchmark 
strategies 

We also computed other statistics in table [6} 



Empirical distribution of performance 




■Uniform 
Optimal 



Performance 



Figure 5: Test 2: Strategy empirical distribution 
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Quantity Formula Value 

Winning percentage 

Relative Optimal Utility 
Relative Optimal Performance 
Utility Sharpe Ratio 
Performance Sharpe Ratio 

VaR 95% Naive Strategy 
VaR 95% Uniform Strategy 
VaR 95% Optimal Strategy 
VaR 90% Naive Strategy 
VaR 90% Uniform Strategy 
VaR 90% Optimal Strategy 



Table 6: Test 2: Other statistics on performance of optimal strategy 

Behavior Analysis In this paragraph, we analyze the behaviour of the strategy as fol- 
lows: first, we plotted in figure [6] the empirical distribution of the number of trades for one 
trading session. Secondly, we plotted trades realizations for three days of the BNPP.PA 
stock for the week starting on 04/19/2010. 

The three following graphs represent three days of market data for which we computed 
the mid-price (lines) with associated trades realizations for the optimal strategy (vertical 
bars). A positive quantity for the vertical bar means a buying operation, while a negative 
quantity means a selling operation. 

Figure [7] shows the trade realizations of the optimal strategy for the day 04/19/2010 on 
the BNPP.PA stock. The interesting feature in this first graph is that we see two buying 
decisions when the price goes down through the 54.5 € barrier, and which corresponds 
roughly to a daily minimum. The following selling decision can be viewed as a failure. On 
the contrary, the two last selling decisions correspond quite precisely to local maxima. 

Figure [8] shows the trade realizations of the optimal strategy for the day 04/22/2010 
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Figure 6: Test 2: Empirical distribution of the number of trades 
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Figure 7: Test 2: Strategy realization on the BNP.PA stock the 04/19/2010. 
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BNP.PA 22/04/2010 




-1000 



Figure 8: Test 2: Strategy realization on the BNP.PA stock the 04/22/2010. 



on the BNPP.PA stock. The interesting feature in this realization is that it looks like a 
U-shaped pattern of liquidation that appears in [lj for an mean- variance optimal strategy 
with a power-law order book resilience. This pattern is very robust, so we expect that it 
appears frequently in the optimal strategy. Moreover, we can notice that the volume of 
each trade in the day is roughly increasing (in absolute value) with the current price. 

Figure [9] shows the trade realizations of the optimal strategy for the day 04/23/2010 
on the BNPP.PA stock. This realization is another illustration of the phenomenon of 
"maxima detection" that appears when executing the optimal strategy. Note that in this 
last realization, the naive strategy was overperforming the optimal strategy, due to an 
unexpected price increase. Despite this, it is satisfactory to see that there are only three 
trades, which is less than on April 19 and 22, 2010, and that the detection of favorable 
price is accurate. 

5.4 Test 3: Sensitivity to Bid/ Ask spread 

In this last section, we are interested in the sensitivity of the results to the bid/ask spread, 
determined here by the two parameters K a and k\>. More precisely, we look at the dominant 
effect between the spread and the multiplicative price impact through the parameter A. 
We proceeded to two tests here: one without bid/ask spread, i.e. K a = = 1 and with 
A = 5.10 -4 as before, and one with a spread of 0.2% and a price impact parameter A = 0. 

Parameters The table [7] shows the parameters of the two tests. We only changed the 
impact and spread parameters and let the others be identical. 
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Figure 9: Test 2: Strategy realization on the BNP.PA stock the 04/23/2010. 



Parameter 


No spread test 


No impact test 


Parameter 


No spread test 


No impact test 


Maturity 


1 Day 


1 Day 


X 


20000 


20000 


A 


5.00E-04 





Yq 


2500 


2500 




0.2 





Po 


51 


51 


7 


0.5 


0.5 


%min 


-20000 


-20000 


Ka 


1 


1.001 


%max 


200000 


200000 


Kb 


1 


0.999 


Umin 








€ 


0.001 


0.001 


Umax 


5000 


5000 


b 


0.01 


0.01 


Pmin 


49 


49 


a 


0.25 


0.25 


Pmax 


53 


53 








n 


30 


30 








rn 


40 


40 








N 


100 


100 








Q 


10 5 


10 5 



Table 7: Test 3: Parameters 



Performance Analysis In tablelSlwe computed several statistics on the results. In figure 



10 we plotted the empirical distribution of performance in the two tests, with the test 2 



distribution (Cf. figure [5]) serving as a reference. In figure 11 we plotted the empirical 
distribution of the number of trades in the two tests, which is particularly helpful for 
interpreting the results. Indeed, we see that a big spread reduces the number of trades 
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of the optimal strategy. This can be interpreted as a phenomenon of clustering, which is 
consistent with the financial interpretation: if the spread is very large, the opportunity to 
buy at low price and sell at high price is significantly reduced, and more risky. Moreover, 
a large spread penalizes strategies that both buy or sell frequently. Then the optimal 
strategy tends to execute bigger quantity in a single trade, supported by the absence of 
price impact. This has the side effect of enlarging the distribution of the optimal strategy 
in this case, since the number of trades tends to decrease. Indeed, we observe that the 
smallest standard deviation is obtained with the uniform strategy, and the largest one 
is obtained with the naive strategy, then, qualitatively speaking, we expect the standard 
deviation of performance to decrease with the number of trades. On the other hand, setting 
the spread to zero does not change the shape of the empirical distribution, and comparing 
table [8] with tables [5] and [6j we see that there is almost no change between zero spread 
and a one-tick spread. This result may be interpreted as follows: first, we notice that 
there is an approximation in the state space, and particularly in the price grid, due to our 
discretization, which is bigger than the scale of one tick (one tick is the price discretization 
unit in the market's limit order book, typically 0.01=6). One finer study would be to set 
the prices grid precisely on the market's prices grid. However, as mentioned before, the 
special feature of our strategy is the path-dependency. Let us consider the typical scale of 
quantities involved in our optimization: we expect the optimal strategy to profit from price 
variation at the scale of 1=6 in our example; if the spread is about 0.1=6, like in our last 
example, and if we usually do about 10 trades on the liquidation period, the effect of the 
spread (10 x 0.1 6= 16) is at the same scale as the price fluctuation. Then, the spread will 
have an important penalizing impact on the optimal trading strategy, and particularly on 
its schedule, i.e. the trading times (the quantities traded are more specifically constrained 
by market impact). On the contrary, if the spread remains small compared to the price 
fluctuations, the optimal trading schedule is not really modified, and the aspect of the 
performance distribution has a similar form. 



Quantity 


No spread test 


No impact test 


No spread vs. T2 


No impact vs. T2 


Mean Utility 


1.00113 


1.00025 


-3.00.10" 5 


-9.08.10" 4 


Mean Performance 


1.00227 


1.00053 


-5.98.10" 5 


-1.80.10" 3 


Standard Deviation 


0.00432 


0.00906 


-9.17.10" 3 


1.078 



Table 8: Test 3: Statistics. In the two last columns "No spread vs. T2" (resp."No impact vs. T2") 
are shown the relative values of "No spread" test (resp. "No impact" test) against the values of test 
2 of the preceding section. 
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Empirical distribution of performance 




-Test 2 
-no impact 



Figure 10: Test 3: Empirical distributions of performance 
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Figure 11: Test 3: Empirical 



Empirical distribution of number of trades (w/o impact) 
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